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Abstract The phenomenological Green's function de- 
veloped in the works of Yang, Rice and Zhang has been 
very successful in understanding many of the anoma- 
lous superconducting properties of the deeply under- 
doped cuprates. It is based on considerations of the 
resonating valence bond spin liquid approximation and 
is designed to describe the underdoped regime of the 
cuprates. Here we emphasize the region of doping, x, 
just below the quantum critical point at which the pseu- 
dogap develops. In addition to Luttinger hole pockets 
centered around the nodal direction, there are electron 
pockets near the antinodes which are connected to the 
hole pockets by gapped bridging contours. We deter- 
mine the contours of nearest approach as would be mea- 
sured in angular resolved photoemission experiments 
and emphasize signatures of the Fermi surface recon- 
struction from the large Fermi contour of Fermi liquid 
theory (which contains 1 + x hole states) to the Lut- 
tinger pocket (which contains x hole states). We find 
that the quasiparticle effective mass renormalization in- 
creases strongly towards the edge of the Luttinger pock- 
ets beyond which it diverges. 
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1 Introduction 

A first qualitative understanding of the superconduct- 
ing properties of the cuprates at optimal and overdop- 
ing is provided by BCS theory with the important mod- 
ification that the gap has d-wave symmetry which re- 
quires nodes and a change in sign on crossing the nodal 
direction, (7r,7r). In the underdoped regime, however, 
anomalous properties are observed which require new 
elements [T2,3,4,5,EllZliaini[J0TlT2^ for their under- 
standing, which fall outside of c?-wave BCS p^lMllTSl 
23 and go beyond effects of inelastic scattering . [TtIITSI 
ITM2U1I^[^ strong coupling, ^511^1^ and anisotropy. [ 
[7711^ Models that might account for some of the un- 
conventional properties have been proposed which in- 
clude the possible role of phase fluctuations, locally 
fluctuating antiferromagnetic order, [3D] and the exis- 
tence of preformed pairs. [3T] These pairs are envisaged 
to form at a higher temperature, T* (the pseudogap en- 
ergy scale), than the superconducting Tc where phase 
coherence is achieved. There are also models of compet- 
ing orders, such as spin or d-density waves. [5^155] Here 
we consider instead the recent model of Yang Rice and 
Zhang (YRZ)pP which is based on a microscopic the- 
ory for a resonating valence bond (RVB) spin liquid. [51] 
YRZ provides an ansatz for the coherent part of the 
charge carrier Green's function. The proposed self en- 
ergy, which describes the formation of a pseudogap on 
the antiferromagnetic Brillouin zone (AFBZ) boundary, 
for doping, x, below a quantum critical point (QCP) at 
a critical doping of x = Xc, has the very desirable prop- 
erty that it is relatively simple and therefore calcula- 
tions of properties are often straightforward. Another 
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important element of the theory, in its final form, are 
the Gutzwiller factors which describe the narrowing of 
the electronic bands as the Mott transition to an in- 
sulating state is approached with reduced doping. An 
additional Gutzwiller factor accounts for a reduction 
in the magnitude of the coherent part of the Green's 
function as the incoherent part grows. 

Many of the superconducting properties of the un- 
derdoped cuprates, which were previously considered 
to be anomalous, have recently been understood 0133 
[551I371I38II51 I1] qualitatively within the model of YRZ 
and this encourages us to explore further its predic- 
tions. The properties considered so far include: Raman 
scattering which shows two gap scales. [2lf35] which 
in the YRZ model are identified as combinations of 
pseudogap and superconducting gap peaks. The large 
reduction in the normalized specific heat jump at Tc 
as X is decreased [7] is also easily understood, [3^ as 
is the penetration depth data[S] indicating that the 
slope of the low temperature, T, linear in T law is 
much less affected by doping than is the magnitude 
of the zero temperature superfluid density. |37j A com- 
parison of Fig. 7(b) of reference [37] with data on Bi- 
2212 for Tc vs the inverse square of the penetration 
depth at zero temperature, which is referred to as the 
Uemura plot, shows good agreement. A critical ele- 
ment of this success is the appearance of the square 
of the Gutzwiller factor, gt{x), in YRZ which goes like 
(j^)^. The data only goes down to transitions tem- 
peratures of approximately one half that of the opti- 
mal critical temperature. It may be that to go closer 
to the end of the superconducting dome in the highly 
underdoped regime, additional effects not included in 
YRZ may need to be considered such as phase fluc- 
tuations. Hat-like structures in the optical self energy 
which emerge in the underdoped regime 38^ and are not 
present at optimum and overdoping are explained as 
is the second energy scale beyond the superconduct- 
ing scale which is seen in the partial optical sum. [3^ 
The checkerboard patternfS^ observed in Fourier trans- 
form scanning tunneling spectroscopy (FT-STS) |40ll4n 
WI\ has been discussed, as have some aspects of the 
angular resolve photo-emission (ARPES) datapi lTUlITT] 
and its relationship to FT-STS 4 . 

An important success of the YRZ model, among 
many others in addition to those just mentioned, is 
that the ARPES results of Kondo et al.[I2] for highly 
underdoped Bi-2201 are easily understood. fBS What is 
measured is the energy of nearest approach to the Fermi 
level as a function of angle, 9, in the Brillouin zone. The 
essential result is that around the nodal direction in the 
non-superconducting state, there exists an arc which is 
ungapped and thus defines a true Fermi surface. At an- 



gles, 9, smaller than a critical angle, 9h, which defines 
the end of the Luttinger pocket, states of zero excita- 
tion energy cease to exist. There is instead an effective 
pseudogap which continues to increase as the antinodal 
direction is approached at small angles. In the super- 
conducting state a superconducting gap is observed to 
form on the ungapped Luttinger arcs and the pseudo- 
gap region remains relatively unaffected by the super- 
conducting transition. This behavior follows directly in 
YRZ theory wherein the pseudogap leads to Fermi sur- 
face reconstruction from a large Fermi liquid (FL) tight 
binding Fermi surface (TBFS) to Luttinger pockets. 
These pockets are centered about the nodal direction 
and are bounded on one side by the antiferromagnetic 
Brillouin zone boundary. However, only the side of the 
Luttinger pocket furthest from the AFBZ boundary is 
strongly weighted and this forms the well known Fermi 
arcs while the backside has a smaller weight and can 
be ignored for many purposes. In fact, it was not until 
very recently that this part of the Luttinger pocket was 
detected in ARPES experiments which previously only 
found the Fermi arc. |43| 

In view of the many successes of YRZ theory, as 
applied to the underdoped cuprates, it is of interest to 
examine more closely a regime which has so far received 
very little attention from both theory and experiment. 
When the pseudogap is large, corresponding to deep un- 
derdoping, one has a Fermi arc around 9 ^ tt/A which 
is followed by a gapped contour of nearest approach en- 
ergy. When the doping is increased towards x = Xc, cor- 
responding to the QCP of the theory, the large TBFS of 
FL theory must be recovered. It is the transition region 
between the hole Fermi pockets, which are characteris- 
tic of the approach to the Mott insulating state, and the 
TBFS of FL theory that we wish to study. This transi- 
tion proceeds through the nucleation of small electron 
pockets bounded by a surface of zero excitation ener- 
gies which appear first in the antinodal direction and 
are additional to the Luttinger hole pockets. There are 
no zero energy states on the bridge connecting these two 
pockets, but rather there is only a contour of nearest ap- 
proach energies on which there is a finite pseudogap. As 
doping is increased further, the gapped contour shrinks 
in size and eventually, both pockets meet, at x — Xc, 
at which point the TBFS is recovered. Note that the 
pockets have two sides, but the YRZ model transfers 
weight between the sides as the pseudogap is reduced, 
creating a continuous transition to zero pseudogap. In 
that limit, the weakly weighted sides define the AFBZ 
boundary. It is the detailed study of the evolution of 
the electronic structure from a large Fermi surface of 
Fermi liquid theory to hole Luttinger pockets that is 
the topic of this paper. We seek here to provide cal- 
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culations that should serve as a guide to experiments. 
Since ARPES provides us, in principle, with a direct 
probe of the Fermi surface evolution with doping, we 
will present our results with this probe in mind. How- 
ever, the Fermi surface reconstruction, which is most 
pronounced in the underdoped region just below the 
QCP will have an effect on other properties as well. In 
particular, we note the possibility of a connection with 
quantum oscillations associated with electron and hole 
pockets. Further, it has been shown to be possible [H] 
to extract from the electronic density of states, N{uj), 
which is measured in scanning tunneling spectroscopy 
(STS) information on the partial contribution to N{uj) 
of each angle, 9, in the Brillouin zone. This type of 
analysis is very similar to ARPES. li-JF] 

Section[2]will describe the relevant theoretical frame- 
work and give the parameters involved in the YRZ cal- 
culations. Section [3] will discuss the electronic structure 
of the pseudogap onset as it should appear in ARPES 
and section U] includes our description of pocket for- 
mation in the context of mass renormalization. Finally, 
section [5] will contain a summary of our conclusions. 

2 Theory 

The YRZ modelfUHl employs, for the coherent piece, a 
Green's function, 

9tW^ 



G{k,uj,x)^ J2 



(1) 



which has been formulated from studies of an effec- 
tive Hamiltonian in an RVB-type, Gutzwiller projected 
mean field t-J model. Entering this propagator are the 
quantities: 



- . a 



a = 



(2) 



The energy dispersion, — —2t{cos k^a + cos ky a) — 
it' cos k^a cos kyU — 2i"(cos 2kxa + cos 2kya) — fJ,p, is the 
third nearest-neighbor tight-binding energy dispersion, 
while ^2 ~ —2t{cos k^a + cos ky a) is the first nearest- 
neighbor term. These energy dispersions contain dop- 
ing dependent coefficients: t{x) — gt{x)to + 3gs{x)Jx/8, 
t'{x) = gt{x)t'Q, and t"{x) = gt{x)tQ, where gt{x) = 
and gsix) = f^i^^y are the energy renormalizing 
Gutzwiller factors for the kinetic and spin terms, re- 
spectively. gt{x) also appears in Eq. ([T]) as a weight- 
ing factor for the coherent part of the Green's function 
which acts to statistically remove or project out doubly 



occupied states. The dispersion, ^fe, uses /Zp as a chem- 
ical potential determined by the Luttinger sum rule in 
the form, 



1 



(2^)2 



d^k. 



G{k,uj=0)>0 



(3) 



which contains the pseudogap in the dispersion of the 
propagator, Eq. (IT|). Thus the chemical potential be- 
comes pseudogap dependent. Values of other parame- 
ters in the dispersion were taken from Ref. 1 to be: 
ioAo = -0.3, t'f^/ta = 0.2, J/to = 1/3, and x = 0.338. 
While here for simplicity we will keep the tight bind- 
ing parameters introduced in the paper of YRZ, they 
could be changed to agree better with the measured 
dispersions in a particular material of interest. 

In the YRZ model, both the superconducting gap, 
Z\sc, and the input pseudogap, Z\pg, have a d-wave k- 
space dependence described by: 



A. 



2 



(cos k^a — cos kytt), 
(cos k^a — cos kytt) , 



(4) 
(5) 



where a is the lattice constant. A'^^{x) and Ap^{x) are 
described by the well known superconducting dome and 
pseudogap line respectively, the latter vanishing at T = 
at a QCP in this model. These are given explicitly as 



A'^^ix) = 0.14io(l - 82.6{x - 0.16)^), 
AlJx)=3to{0.2-x). 



(6) 
(7) 



We have used the well known empirical expression of 
Eqn. (|6]) for the Tc dome with the assumption that 
twice the gap to Tc ratio is 6. Recently, Schachinger 
and Carbotte [15] have solved a BCS equation for gap 
and critical temperature in the YRZ model and have 
found that the opening of a pseudogap has the effect 
of increasing this gap to Tc ratio as the doping is re- 
duced. This implies that the superconducting gap re- 
mains large in underdoped cuprates, consistent with 
their small coherence length. 

From the YRZ Green's function of Eq. ((Ij , one can 
extract the spectral function, A(k,io), given by 



A{k,cu) = J2 gtW^[iun'Siu;-E^) + iv'^fdicu + E:)], 

(8) 

where 

(9) 



E-^^{E-y + Ai 









[1(1 


-f)] 





(10) 

(11) 
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Here, E" represents the pair dispersion and (m")^ and 
(w")^ are standard Bogoliubov quasiparticle amplitudes 
for their respective a energy branches. 

3 ARPES 

Central to the YRZ model is the existence of electron 
and hole pockets. These regions of occupied or unoccu- 
pied states overlap in the YRZ model with weighting 
factors of which enter much the same as in BCS 
theory (Bogoliubov and hole-particle weightings) , 
but here, Age is replaced with Apg and replaces the 
quasiparticle energy as in Eq. ([2]). While and are 
centered about the Fermi level, the W'^ and W~ are 
centered about the zeros of a different energy, (Cfc + 
whose zeros acts as the origin of the pseudogap. This 
origin corresponds to the TBFS, given by = 0, be- 
ing shifted by an effective particle kinetic energv. |T7) 
which near half filling corresponds to the AFBZ. In the 
Apg — > limit we therefore establish two boundary en- 
ergies in this system, that of the TBFS, and that of the 
AFBZ, wherein have values dependent upon the 
sign of ^fe + Cfc- In this limit the TBFS has a weight 
(W^) of 1 and the AFBZ a weight of zero; thus it is 
the weighting factors which provide the agency whereby 
the AFBZ drops out of the quasiparticle energy in this 
limit, as we know it must. 

A fundamental aspect of the work presented here is 
the determination of nearest approach contours and en- 
ergies. Wc define these quantities by marking local min- 
ima of the and E/^ energies in the Brillouin zone, 
which we will refer to as E^^, and that occur at loca- 
tions kna- Often it is convenient to display the energies 
as a function of angle, 9, as measured from the corner 
point at (tt, tt). This allows us to compare to a simple 
d-wave superconducting gap such as Asc{0) « cos(26'). 
The angle and magnitude of k defines the momentum in 
the first Brillouin zone. For some angles there are mul- 
tiple local minima due to the hole and electron pockets; 
but also, there are local minima which bridge the pock- 
ets. These connections, or 'bridges', do not appear in 
other contours such as constant energy surfaces nor in 
the Luttinger surfaces. In fact, it is for angles along 
these bridges that the pseudogap has a finite, effective 
pseudogap value. Such a set of contours is shown in 

Fig.m 

On the E^^ energies, which arc > for all an- 
gles, there exists a zero energy Fermi surface for angles 
9 < 9e, which surrounds a region of filled electron states 
(the 'electron pocket') near (0,7r) and (7r,0) and is shown 
along with the nearest approach contours in Fig. For 
the E~^ energies, which are < for all angles, there ex- 
ists a zero energy Fermi surface for angles 9 > 9h, which 




Fig. 1 (a) Example of heavily weighted nearest approach con- 
tours for X = 0.12 and x = 0.18. The influence of the 
band can be seen in the x=0.18 case near {tt,0) and (0,7r) as 
the appearance of electron pockets, (b) Full nearest approach 
contours for E+ and in the lower right portion as well 
as the zeroes of in the upper left portion, (c) Nearest ap- 
proach pockets for E+ and E~ (here equivalent to Luttinger 
surfaces) for a single doping of i'=0.18 for varied input ^pg. 
The pseudogap, zipg = O.OQto , corresponds to YRZ fitted pa- 
rameters for X = 0.18 as described in the text. 

surrounds the region of unoccupied states (the 'hole' or 
Fermi pocket) along the AFBZ boundary. Important 
evidence of a Fermi pocket has been obtained exper- 
imentally by Meng et al. .43, on Bi2Sr2_a;La2;Cu06+5, 
which revealed not only the presence of a Fermi pocket 
but also its coexistence with an additional gapped sur- 
face extending the Fermi arc beyond the corners of the 
Fermi hole pocket. This experimental observation is an 
expected feature of the YRZ model as the manifestation 
of these bridge contours in ARPES data. 

We depict this behavior in Fig. [IJa) where we have 
shown the most strongly weighted nearest approach 
contour with energy < 0. For angles near 45° we display 
the back of the Fermi pocket. For angles 9 < 9h, which 
marks the corner of the hole pocket, we display the con- 
tour of nearest approach for the E~ band, which has 
negative energies, as shown in Fig. [5] In this region, the 
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Fig. 2 i?+ and E energies vs along their individual con- 
tours of nearest approach. Shown for x = 0.19, 0.195, 0.199 
and 0.20 

nearest approach energy is a local minima only, caused 
in the YRZ model by proximity to a zero in (^^ + ^2) 
which corresponds to a zero valued minima in This 
choice of contour is maintained until 9 < 9^ ai which 
point the zero energy contour of the E'^ band (shown 
in Fig. becomes the more strongly weighted contour 
of nearest approach. One will notice that the contour 
may have kinks when crossing either of 9^ or 9h ■ 

The choice of contours in Fig. [ija) is justified on 
the criteria that the energies are: 1) a local minima in 
the absolute value of energy, 2) of negative (or zero) 
energy, 3) have the largest weighting or W^. The 
full set of contours, for a;=0.18, as well as the (^^ + 

= surface are shown in Fig. [T](b). It is precisely 
this effective kinetic energy contour given by = 
which forms the bridged regions between the pockets 
and therefore dominates in the well underdoped, when 
the hole pockets have become very small, and behaves 
as an effective Fermi surface in that case. 

It can be seen in Fig.[ljc) that changing the magni- 
tude of the pseudogap directly modifies the size of the 
Fermi pocket. However, there are two 'hard' barriers in 
the structure of the hole pocket. The first is the AFBZ 
boundary. The second is, in fact, the normal state Fermi 
liquid TBFS. As the pseudogap value decreases, the 
pocket attempts to stretch out to merge with the TBFS. 
However, for small pseudogap, the hole pocket is re- 
stricted to the bounded region created by the AFBZ 
and the TBFS. Further expansion in area is not al- 
lowed, as the Luttinger sum rule dictates, in essence, a 
conservation of the ratio of hole/electron pocket areas. 
Therefore, beyond some pseudogap value, for a given 
doping, electron pockets must emerge in the corner re- 
gions (which are also bound by the TBFS and AFBZ) 
which extend into the upper antiferromagnetic Brillouin 
zone. The area of these electron pockets balances the 
increased area of the hole pocket, finally merging at the 
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Fig. 3 Weighting factors along nearest approach contours. 
For the cases with no superconductivity, the weightings are 
iy+ and W~ . The inclusion of superconductivity adds an 
additional factor of which represents the occupied state 

weighting of the Bogoliubov quasiparticles. 

intersection of the AFBZ and TBFS for Z\pg -> 0. We 
note that the heavily weighted hole and electron pock- 
ets for Z\pg = 0.06 evolve into the sohd red contours for 
Z\pg = 0.01 with the lightly weighted sides collapsing 
onto the antiferromagnetic Brillouin zone. In the limit 
of Z\pg = this part has zero weight. This shows how 
the energy which plays a role in YRZ, drops out 
of consideration in the dispersion curves in the limit of 
the Fermi liquid, where only plays a role. 

Examples of the weighting factors along nearest ap- 
proach contours are shown in Fig [31 The top frame is 
for x=0.12, where there are no electron pockets, and the 
bottom frame is for x=0.18 which lies between Xonset 
and Xc and thus has electron pockets. For the x=0.12 
case, is single valued, while for x=0.18, the pres- 
ence of an electron pocket splits the contour at small 
angles. The dashed black line, 9 < 9^, represents the 
side of the electron pocket nearest to the AFBZ. The 
weightings are double valued for both dopings and 
are shown in red, where the dashed curve is now rep- 
resentative of the side of the hole pocket nearest to the 
AFBZ. It is only at 9 = 7r/4 that the sohd red fine 
reaches its Fermi liquid value of exactly one. This oc- 
curs even though there exists states of zero excitation 
energy on the entire hole Luttinger Fermi surface just 
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as in an ordinary Fermi liquid. Here, < 1 represent 
a departure from FL theory. This departure increases 
particularly quickly as the end of the hole pocket is ap- 
proached where W~{6) has already dropped below 1/2 
in value for both dopings shown. To understand how 
this comes about, we note that on the hole Luttinger 
surface Ef, = and that the weighting factor on 
this contour simplifies to 



1 - 
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VSna 



VSna / 



(12) 



On the front part of the hole Luttinger pocket fc„o is 
away from the antiferromagnetic zone boundary, which 
is defined by = 0, but on the opposing side, we are 
close to the AFBZ boundary and so we can expect 
to be small. This means that on the front side of 
Eq. [T2]is of order 1, but is reduced due to the presence 
of the gap, Z\pg(fc„a), which is squared in the denomina- 
tor. Right at 6 = it/ 4 the value of Wj^ reaches exactly 
one as shown in Fig. [3] (solid red curve) and W,^ has 
dropped to a little less than 1/2 its value for 9 = Oh, 
marking the end of the hole pocket. On the other hand, 
the weight on the backside (dashed red curve) rapidly 
drops from its maximum a.t 6 — 9h as 9 is increased 
towards the nodal direction. As the nodal direction is 
approached, the contour fc„Q shifts closer and closer 
to the AFBZ boundary and ^^{kna) becomes smaller 
and smaller and so of Eq. [T^] rapidly drops to a 

value near zero. On the gapped contour beyond 9h for 
a; = 0.12 (top frame) and in the bridging region between 
9h and 9^. for a; = 0.18 (bottom frame) the last equal- 
ity in Eq. [T2]no longer holds, but numerical evaluation 
shows that Wy, continues to drop, reaching a value 
« 0.3 in both cases. These behaviors are quite distinct 
from the Fermi liquid case of a large Fermi surface for 
which the weight would remain 1 at all angles. 

The behavior of associated with the nearest 

approach contour of the E'^ (solid black line) is most 
interesting for x = 0.18 (bottom frame). In this case, 
for 9 between and 9e we have two more zero energy 
contours; one with large weight (solid line) and one with 
smaller weight (dashed line) . The weights on these sur- 
faces, however, remain intermediate, well above and 
well below 1. It is in this respect that these electron 
pockets maintain continuous electronic states which sup- 
port closed orbits within a single band unlike the hole 
pocket, which has no continuous state in the nodal di- 
rection near the AFBZ. 

There are complications involved in the inclusion 
of superconductivity. In this case we must also include 





(a) 









Fig. 4 One quarter of Brillouin zone displaying: Top Row: 
and E'^ in (a) and (b) respectively, for the a;=0.18 case 

where the color scale has units of to. Bottom Row: and 
, (c) and (d) respectively, also for the x=0.18 case. 




20 30 
e (deg) 

Fig. 5 Top frame: Ap^J are the nearest approach energies of 
the most strongly weighted of E~ or i?+ bands, for various 
dopings. For large pseudogap, the system is dominated by 
the E^ band. At x=0.17 the i?+ band provides an i? = 
Fermi surface at small angles with a strong weighting. Bottom 
Frame: Same as above but for i?^. This gives an additional 
contribution from Aac- Note for comparison to ARPES plots, 
we have taken the absolute value of these energies. 
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the Bogoliubov u| weightings, which take the form of 
Eq. pT|) We show the W^^^y^^ (sohd green curve in 
Top frame of Fig. [3]) and W'^v^^" (sohd blue curve in 
bottom frame of Fig. [3]) , along the contours shown in 
Fig. [Tfa). On the surfaces defining the hole and elec- 
tron pockets, and E'^ are respectively zero, so that 
on both the front and back side of these surfaces, the 
Bogoliubov is precisely 1/2. This is shown as the 
green and blue sohd lines in Fig. [TJa). For 9 < Oh and 
in the interval Oh to Of. for x = 0.12 and x = 0.18 re- 
spectively, the reduction in Wv from W alone is less 
than 1/2 because on these contours E^^^ is negative 
and finite. Hence, v~'^ in Eq. ([TT|) is greater than 1/2. 
For the specific case of x = 0.12 (top frame) the solid 
red and solid green contours have practically merged at 
= while in the lower frame, the solid blue contour 
increases towards the solid red until = Oe at which 
point we change from the bridging contour to the new 
electron surface associated with the electron pockets. 

For completeness, we illustrate the E^ energies and 
their weighting factors in the full Brillouin zone in Fig. 01 
The electron pockets are shown in white near the antin- 
odes of Fig. IH^a), while the hole pocket is located in 
the nodal region and is shown in orange in Fig. lU^b). 
The spectral functions of the two dispersions are over- 
laid with their respective weighting factors shown in 
Fig. |3Kc)and (d). It should be clear that the weighting 
of the E~ bridges should be reduced as compared to 
the Fermi pocket itself. 

Recent experimental determination of gap values in 
Bi-2201 using ARPES [12] have been qualitatively de- 
scribed by the YRZ model for the well underdoped 
case where the pseudogap dominates the system |35j. 
Here we present similar calculations over a full dop- 
ing range and demonstrate the additional features cre- 
ated by electron pockets in the small pseudogap case. 
Fig. [SJa) shows this effect in the absence of supercon- 
ductivity. For X = 0.20 we would obtain a full zero 
energy, gapless Fermi surface. For small dopings, say 
X = 0.05, we have a strong pseudogap which follows 
a modified cos(20) function, terminating at Oh ~ 35°. 
This qualitative structure is maintained until the dop- 
ing Xonset where the electron pocket has now formed a 
zero energy TBFS with strong weighting at small an- 
gles, < Of., similar to the energies shown in Fig.[5J 
Therefore, for dopings where Xonset < x < Xc, the curve 
appears as a modified cosine function for angles > Oe, 
and drops to zero for < Oe- This leads to the 'russian 
doll' type arrangement of dashed curves covering the 
four largest values of x considered in the figure. This 
feature is still present when including a superconduct- 
ing gap, as is shown in Fig. [HKb). The extreme cases 
of Age = or Apg = produce what can be approx- 
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Fig. 6 Progression of hole and electron pockets, defined by 
their angles of onset, de and 6^ as a function of doping con- 
verging to e* ^ 10° a.t X = Xc = 0.20. 

imated as cosine functions terminating at Oh and 45° 
respectively. Again, there is a narrow region of doping 
between Xonset and Xc wherein the behaviors of these 
curves are modified. 

For < Oe and > Oh, the behavior is that of the 
superconducting gap. This corresponds to the contour 
of the normal state Fermi surface where the supercon- 
ducting gap is dominant. For Oe < < Oh the behav- 
ior is that of the finite valued, but negative, E~ band 
energy which is interpreted as an effective pseudogap 
value. This corresponds to the pseudogapped out re- 
gion between the hole and electron pockets which fol- 
lows along the nearest approach bridges where there is 
no Fermi surface in the normal pseudogapped state (see 

Fig.mb)). 

Of primary interest for this work is the behavior of 
the pocket edges. Fig. [H] depicts the calculated Oe and Oh 
critical angles across all dopings and their convergence 
to the value 6** at a; = 0.20. This figure is separated 
into two regions for < 0* and > 0* for the elec- 
tron and hole pocket angles respectively. The angle Oe, 
which is always less than 0* , maintains a zero value cor- 
responding to the absence of an electron pocket until 
X = Xonset abovc which Oe is finite, while Oh, which is al- 
ways greater than 0* , has finite value at all relevant dop- 
ings. The angles Oh and Oe converge at a; = Xc = 0.20 
to a value of 6' « 10°, which corresponds to the inter- 
section of the TBFS and the AFBZ. The doping de- 
pendence of Oh is nearly linear away from Xonset , show- 
ing two distinct regions, x < x onset and x > x onset, 
which differ in slope as illustrated by the black dashed 
lines which intersect at Xonset ■ The modification to the 
slope coincides with the onset of electron pockets. Ex- 
perimental verification of the modification of this slope 
would be evidence of the presence of shallow electron 
pockets forming as dictated by the Luttinger sum rule. 
We hope our work will stimulate such experiments. 
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4 Effective Mass 

We wish to discuss the effective electronic mass on the 
Luttinger surfaces in the pseudogapped state. In the 
normal pseudogap state, one can reformulate Eq. [T] as 
it was presented in Ref. pQ by using the pseudogap self 
energy which is given by: Z'pg(fc,a;) = zipg(fc)/[a; + ^2], 
which depends on momentum fe as well as on energy, 
w. The inclusion of a pseudogap self energy acts to 
directly renormalize the energy of the quasiparticles. 
These renormalized quasiparticles will behave as if hav- 
ing a new renormalized effective mass, which we now 
calculate. To begin this, for a; — > 0, the Green's func- 
tion has a simple pole at 



-\ -1 



(13) 



We evaluate Eq. on the Luttinger surface which 
defines a contour in momentum space k = [kp,6] where 
kp is the magnitude of k at fixed angle 6. In the case 
of multiple kp for a given 9, we are interested in only 
the strongly weighted portion, as displayed in Fig. |31 
However, on this contour, Ej^ of Eq. ([2]) is zero and 
therefore 



(14) 



For k near kp at fixed 9 we use k ~ kp + Sk. As a result, 
w in Eq. (jl3p can be rewritten as 




Z\2 

Pg 



(fc) 



0\2 



h—kf 



5k, 



(15) 



where use was made of Eq. 

In the underlying Fermi liquid state, where Z\pg(fc) = 
0, the corresponding quasiparticle pole with which we 
wish to compare would be at = ^P- 5k. The 

k=kF,0 

ratio of these poles gives the ratio of renormalized to 
unrenormalized electronic mass coming from the pseu- 
dogap formation as: 



m;(Z\pg = 0) 



Alg{kF,0) 



_d_ 

dk 



dk 



k=kF.S 



(16) 



Here we have made the approximation that the underly- 
ing Fermi liquid band structure varies only slowly with 
k around the Fermi energy and have evaluated at 



the pseudogap nearest approach contour rather than 
the Fermi liquid contour. This is expected to be a good 
approximation as the two contours do not differ much 
on the strongly weighted side. Fundamentally, we are 
restricted from evaluating this quantity in this simple 
manner on the opposing lightly weighted side of the 
pocket as the regions in k-space differ completely for 
the two cases, with and without pseudogap, causing the 
approximations here to break down. Results for the ef- 
fective mass ratio of Eq. (jl6p are presented in the lower 
frame of Fig. [7] for three values of doping. Before de- 
scribing these in detail, we note two limits of Eq. ([T6|. 
First, in the nodal direction, Z\pg(A;) — and all en- 
ergies reduce to the Fermi liquid unrenormalized band 
structure, in which case mg(Z\pg)/mg(Z\pg = 0) = 1. 
Also, when the hole pocket ends dX 9 = 9h the deriva- 
tive of on the Luttinger contour with respect to |fc| 
at fixed 9 vanishes, ie. -^E~{k, 9h) — for fixed 9 — 9}i. 
This leads to the relation 



(17) 



from which we can conclude, making use of Eq. ([T4l) . 
that the quantity in the denominator of Eq. (jl6p van- 
ishes and so the effective mass renormalization is infi- 
nite. The quasiparticles at the termination of the Lut- 
tinger hole pocket have become infinitely sluggish to 
continued progression m.9. To treat the case when elec- 
tron pockets are present is a straightforward general- 
ization of the above to their respective contours. 

In Fig. [7] we summarize our numerical results for 
three values of doping x = 0.05 (heavily underdoped) 
X = 0.12 and x = 0.18 (which is near our chosen quan- 
tum critical point ai x = Xc = 0.2). In the top frame we 
show the slopes of the electron dispersion curves across 
the Luttinger contour at w = 0. The Fermi liquid case 
corresponds to the dashed curves and the renormalized 
YRZ pseudogap case corresponds to the solid curves. 
The shaded region helps to see at a glance the difference 
between these two cases. As we have already discussed, 
at 9 = 7r/4 (nodal direction) there is no slope renor- 
malization. Away from this point, however, the mass 
renormalization increases and at the end point of the 
Luttinger pocket it diverges as we have already demon- 
strated on the basis of Eqns. and (IT71) . For any 
doping exhibiting electron pockets {x onset < x < Xc) 
there is an additional region about the antinodal di- 
rection which starts at 9 ^ 9e and ends at 6' = 0. At 
9 = 9e the mass is again infinite, as we have zero slope 
for the quasiparticle dispersions after which it gradually 
decreases and has its minimum of about 1.7 at 9 ~ 0, 
in the case shown (x = 0.18), normalized to the Fermi 
liquid case. The value of m*{9 = 0) will shift contin- 
uously to 1 as the pseudogap is reduced in magnitude 
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Fig. 7 a) The slope of crossing k = kp as a function of 

0. The solid curves have pseudogap, while the dashed curves 
are the slopes across the underlying TBFS of ^(fc). Here a is 
determined as in Fig. 3Ja), as the energy branch that main- 
tains coherent quasiparticles. b) The ratio of slopes in (a), 

ii[^;fe„J^Ps)]/i:[^£„J'4pg = 0)1, which is also given pre- 
cisely by the effective mass ratio of Eqn. 1161 

or X approaches Xc- On the other side of x onset, the 
physics of the approach to the Mott insulating state at 
half filling proceeds with a loss of electron states which 
have a zero excitation gap. These reside on the signifi- 
cantly weighted Luttinger arc which decreases in length 
as doping, x, is decreased. In addition, the quasiparti- 
cle effective mass {m*) increases above its Fermi liquid 
value as the end of the arc is approached where m* 
diverges. 

In all cases, however, there remain coherent quasi- 
particles, granted fewer and fewer, as the Mott tran- 
sition is approached, as is also observed in the work 
of Lee et. al. [35] These authors find a coherent opti- 
cal response which remains down to the lowest dopings 
considered, well past the end of the superconducting 
dome in the underdoped cuprates. A second observa- 
tion long considered anomalous, which finds a natural 
but very robust explanation in the YRZ model, is the 
work of Zhou et. al.[in] who found that the slope of the 
dispersion curves in the nodal direction in LSCO are 
only very weakly dependent upon doping. We see here 
that this slope remains unaffected by the formation of 
a pseudogap and hence it is dependent only on the un- 
derlying band structure. The energy scale involved in 



this case is the bandwidth scale or, more precisely, the 
effective first nearest neighbor hopping scale (to)- Mod- 
ification of doping only involves small changes in the 
position of the Fermi momentum with respect to its 
value at half filling, and this will not change the cor- 
responding nodal direction Fermi velocity much, as is 
observed. There also exists evidence [^ [TU1[TT| that the 
length of the ungapped arcs about the nodal direction 
increases in a roughly linear fashion with temperature. 
This could be due to a reduction in pseudogap am- 
plitude with increasing temperature, T, which can be 
justified with the inclusion of scattering. It could also 
be modeled with an input pseudogap which opens more 
slowly with angle out of the nodal direction as one raises 
the temperature. 

Although not rigorously explored here, we take note 
of the consistency of the antinodal effective mass results 
along the small electron pocket in the x=0.18 case. The 
value of m* w 1.7 at = agrees well with the results 
of de Haas- van Alphen (dHvA) and Shubnikov-de Haas 
(SdH) quantum oscillation experiments. These exper- 
iments have reported values of m-^jj^^ — 1.76 ± 0.07 
and Tigjjj = 1.9 ± 0.1 [5U1ICT] . consistent with electron 
pockets [52] ■ Further, the area of the electron pock- 
ets is restricted in YRZ by the underlying TBFS, and 
AFBZ to a value of at most 2 — 3% of the Brillouin 
zone, also in reasonable agreement with experiment. Of 
course, the weightings of Fig. [3] results in a reduction 
of coherent states on the hole pocket near AFBZ which 
could suppress the signal of oscillation from the hole 
pocket. By contrast, the weighting along the electron 
pocket always remains significant in magnitude. There 
may be conditions under which the hole pocket can sus- 
tain coherent oscillations, consistent with a hole area 
{Ah) of percent given by the Luttinger sum rule to be 
Ah — f where subscripts e represents the emerging 
electron pocket. Otherwise, it may be that, with loss of 
fully weighted pockets, oscillations in YRZ should re- 
vert to a Fermi arc picture where it has been shown 
that quantum oscillations could still arise but with the 
frequency of oscillation related to the arc length, rather 
than to the area enclosed by the pockets. |5 3 i Clearly, 
the full calculation of quantum oscillations in the YRZ 
model would produce interesting results. 



5 Conclusions 

We have studied how the electronic structure in the 
cuprates evolves due to the opening of a pseudogap in 
the model of Yang, Rice and Zhang. The pseudogap 
is associated with the antiferromagnetic Brillouin zone 
boundary and has its microscopic origin in a RVB spin 
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liquid. While we considered the entire underdoped re- 
gion of the phase diagram, an emphasis has been placed 
on the region near the quantum critical point at doping 
X = Xc below which the pseudogap becomes finite. We 
also emphasize how the large Fermi surface of Fermi 
liquid theory progressively transforms into a small Lut- 
tinger hole pocket centered about the nodal direction, 
a characteristic feature of the loss of metallicity on the 
approach to a Mott insulating state at half filling. Just 
below Xc (QCP), important modifications occur at the 
point of crossing between the TBFS and the AFBZ, 
where there is no longer a zero energy contour. This 
breaks the large Fermi surface into two pieces separated 
by a bridging region. Both pieces are able to sustain 
zero energy excitations, but the bridging region, be- 
tween a large hole pocket, in the nodal direction, and 
the small electron pocket near the antinodal direction, 
is gapped. Each pocket has two sides, one with a weight- 
ing factor near one and the other with a smaller weight- 
ing closer to zero. The strongly weighted piece is close 
to the large TBFS while the weakly weighted piece is 
close to the AFBZ line associated with half filling which 
corresponds to a Mott insulating state. As the doping 
is further reduced, the hole pocket shrinks in size as 
does the electron pocket which eventually disappears 
entirely leaving only a hole pocket (in each quadrant). 
When it exists, the electron pocket is limited to angles 
in the interval from 9 = to 9 = 9e m the antinodal re- 
gion. A similar situation holds for the hole pocket which 
extends from 6h to 7r/4. 

The behavior of the effective pseudogap as a func- 
tion of angle along the bridging contour is particularly 
interesting. Its value starts from zero aX 9 = 9h and 
grows until 9 — 9^ s.t which point it shows a discon- 
tinuous drop to zero as we move to the electron pocket 
piece of the Fermi surface which is not continuously 
connected to the bridging contour. When superconduc- 
tivity is present there is a superconducting gap on the 
ungapped parts of the normal state Fermi contour, but 
in the region between 9e and 9^, there remains an en- 
hancement, over its superconducting state value, of the 
effective gap measured in ARPES. Of course when we 
approach half filling the electron pocket ceases to ex- 
ist and 9e = 0. Wc find that 9^ is non-zero only above 
X — Xonset ~ 0.167, a value which depends somewhat 
on the assumptions made for the angular dependence 
of the input pseudogap that enters the YRZ model, 
as well as details of the band structure parameters of 
the model. These authors also use the lowest da;2_j,2 
harmonic, Eq. [5l for the momentum variation of the 
pseudogap in the Brillouin zone. If this assumption is 
relaxed and this harmonic is taken to power n then vari- 
ations in the exponent, n, can change the shape of the 



Luttinger contours as well as the onset doping, Xonset- 
These are, however, secondary considerations and the 
results we have just described remain almost unmodi- 
fied for large variations in the input pseudogap angular 
profile; for example, a fc independent pseudogap. Thus 
the physics described here is robust. 

On the Luttinger part of the nearest approach con- 
tours, the effective pseudogap is zero. Nevertheless, the 
energies in this region are modified from those associ- 
ated with the underlying TBFS. A measure of these 
changes is the change in effective mass at fixed angle, 
9, that is brought about by the presence of the input 
pseudogap. The effective mass is unchanged from its 
tight binding value only in the nodal direction. As we 
move towards the end of the Luttinger hole pocket the 
effective mass increases and has a divergence right at 
the corner where we go from linear in momentum, away 
from the Fermi momentum, to a gapped behavior for 
the renormalized dispersion curves at fixed angle, 9. 
A similar effect applies when the doping is sufficiently 
close to the QCP at a; = a;c that an electron pocket 
also exists. In this case, the renormalized mass is larger 
than the bare mass everywhere and diverges a,t 9 — 9e- 

The picture that emerges for low doping, near half 
filling, is an approach to the Mott transition where the 
number of ungapped quasiparticles available becomes 
greatly reduced, residing only on the Luttinger sur- 
face around the nodal direction. In addition, as the 
size of the Luttinger hole pocket shrinks the effective 
mass of the few remaining quasiparticles on the heavily 
weighted Fermi arc retains its tight binding value only 
aX 9 — Tr/4: but increases as the end of the Luttinger 
pocket is approached where there is a divergence. 
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